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Abstract 

Ultra high energy (UHE) particles of cosmic origin impact the lunar regolith and produce radio 
signals through Askaryan effect, signals that can be detected by Earth based radio telescopes. We 
calculate the expected sensitivity for observation of such events at the Giant Metrewave Radio 
Telescope (GMRT), both for UHE cosmic rays (CR) and UHE neutrino interactions. We find that 
for 30 days of observation time a significant number of detectable events is expected above 10 20 
eV for UHECR or neutrino fluxes close to the current limits. Null detection over a period of 30 
days will lower the experimental bounds by a magnitude competitive to both present and future 
experiments at the very highest energies. 
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I. INTRODUCTION 



The study of cosmic rays with energies near and above the Greisen-Zatsepin-Kuzmin 
(GZK) cut-off ~ 5 x fO 19 eV is of great interest in particle astrophysics. Above the 

GZK energy, proton interactions p + 7cmb ~^ Ntc on Cosmic Microwave Background (CMB) 
photons are possible. The universe in effect becomes opaque over Mpc scales. To check 
the existence of this cut-off is therefore an important motivation for cosmic ray experiments 
like AUGER [3| and HIRES 4|. Complementing the cosmic ray experiments, a number of 
experiments perform direct searches for ultra high energy (UHE) neutrinos from the GZK 
pion decays (tt + — > /U + z/^ — > e + v e v^v^). The search for neutrinos above the GZK energy is 
also motivated in the context of Grand Unification Theories (GUT) [5|. Relic 'X'-particles 
with ultra high mass Mx ~ 10 22 eV originating from the GUT phase-transition in the 
early universe could be decaying in the present epoch. Thereby they would produce UHE 
neutrinos. As a final example of a theoretically very appealing possibility, is the Z-burst 
process [&]. In this scenario, the highest energy cosmic rays are produced following resonant 
interactions [yv — > Z°) of UHE neutrinos on the relic neutrino background. To explain the 
full cosmic ray spectrum, this would require a very high flux of neutrinos in the E v ~ 10 22 
regime to be present. Limits constraining the possible UHE ! neutrino flux from these 
processes exist, with the most stringent one to date being presented by the ANITA-lite 
experiment [7(. 

Among the various methods for detecting ultra high energy (UHE) particles, a promising 
technique utilizes the Askaryan effect S]. This effect causes electromagnetic cascades, in- 
duced by the UHE particle interactions in dense media, to develop a negative charge excess. 
The charge excess will radiate coherently at radio wavelengths, thereby producing strong 
coherent pulses of Cerenkov radiation. In a series of accelerator experiments, the Askaryan 
mechanism has been confirmed to work in different media. Most important for our purposes, 
it was tested in Silica sand specifically to mimic the conditions of the lunar material 9j. It 
was first proposed by Askaryan Q], later by Dagkesamanskii and Zheleznyk 10], that the 
radio transparent lunar regolith (the 10 - 20 m deep surface layer of the Moon, consisting 
mainly of fractured rock) would be an ideal target for studying such UHE particles. An 
additional advantage is the absence of an atmosphere which implies that electromagnetic 
showers in the lunar regolith are caused by the primary cosmic particles. 
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A number of attempts at observing cosmic ray induced radio waves from the Moon have 



since been carried out using ground based radio telescopes. The first experiment [11] used the 
64 m diameter Parkes radio telescope in Australia to make coincidence measurements from 
two polarization channels. A 500-MHz band centered at 1.425 GHz was used after making 
appropriate corrections for the ionospheric delay between two sub-bands. The second was the 
Goldstone Lunar Ultra-high energy neutrino Experiment (GLUE) [3| that used two large 
telescopes of the JPL/NASA Deep Space Network in Goldstone (CA). A third experiment 
of this kind took place at Kalyazin Radio Astronomical Observatory, using again a single 



64 m radio telescope 13]. None of these experiments however, resulted in any detection 



of signals from UHE particles. Under the name NuMoon, there is one currently on-going 



search using the Westerbrok Synthesis Radio Telescope (WSRT) in the Netherlands [ijj]. 
For future experiments with increased sensitivity there are proposals 15|, [l6j, ll7] to study 



this effect with upcoming radio telescopes like the Low Frequency Array (LOFAR) and the 



Square Kilometer Array 



more radio antennas 



18 



SKA), or to use lunar orbiting artificial satellites carrying one or 



19]. 



In this paper we study in detail the possibility of detection of UHECR and neutrino radio 



waves from the Moon with the Giant Metrewave Radio Telescope (GMRT) facility 20j. The 
GMRT comprises 30 fully steerable dishes spread over distances up to 25 km, each with a 
diameter of 45 m. Half of these antennas are distributed randomly over about 1 km 2 in 
a compact array while the rest are spread out in an approximate "Y" configuration. The 
GMRT currently operates in eight frequency bands around 150, 235, 325, 610 and 1000-1450 
MHz. The shortest baseline is 100 m while the longest one is 26 km. The telescope has 
an effective area in total of 30000 m 2 for frequencies up to 327 MHz and 18000 m 2 at the 
higher frequencies. The RMS sensitivity of the GMRT, using a 2 s integration time and a 
bandwidth of 32 MHz is 0.3 mJy and 0.03 mJy at 327 MHz and 1.4 GHz respectively. With 
such high sensitivity, we expect a substantial improvement from previous experiments in the 
measurement of the Askaryan radio waves from the Moon. 

II. RADIO WAVES FROM CASCADES IN THE LUNAR REGOLITH 

High energy charged particles from UHE cosmic ray (CR) or UHE neutrino interactions 
in the lunar regolith will initiate a cascading shower with total energy E s and typical length 
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scale [17] 

in units of radiation length. For the lunar regolith with radiation length Xo = 22.1 g/cm 2 
and density p = 1.7 g/cm~ 3 , the shower length for a particle with energy E = 10 20 eV 
is therefore 1.7 m. Since the particles in the cascade travel at speeds very close to the 
speed of light the duration of the shower is L/c ~ 5.67ns. The radio waves of about 6 ns 
duration which are emitted from the lunar regolith are dispersed as they propagate through 



the atmosphere of the Earth. The electrons in the ionosphere cause a time delay 16] of 



i, H = 1.34 x ID"' (^) (^) 2 (2) 

seconds, where N e = 10 17 m -2 is the typical night-time column density of electrons. The 
dispersion of a radio signal of frequency v and bandwidth Av at zenith is therefore 

Au \ /200MHz x 3 



At d = 134ns x . (3) 

V 40 MHz / V v J K ' 

One must thus look for Cerenkov radio pulses of about 100 ns duration from the Moon and 
use coherent de-dispersion to recover the broadband structure of the signals. 

The intensity of radio waves on Earth from a shower in the lunar re goli th with energy E s 



emitting Cerenkov emission has been parameterized from simulations 



13, 



211 ] . At an angle 



9 to the shower axis, for radiation frequency v and bandwidth Au, it is given by 



E* \ 2 ( 1 m\ 2 / v \ z Av 



n*,M = 5-3 x 10- /(*,») x^-J [ Ml + {vM1M] ) iooMiiJyW 

where R is the distance between the emission point on the Moon's surface to the telescope 
and vq = 2.5 GHz. Furthermore there is an angular dependence given by 

sin# x 



f{E.,0)= — - exp(-2Z 2 ) (5) 
\sm V C J 

with Z = ^(cos^c - cos^) and k(E s ) = (z//GHz) 2 (70.4 + 3.4 \n(E s /TeV)). The angular 
distribution of the radiation pattern is illustrated in Fig. [21 In this figure we also show the 
commonly used Gaussian approximation where the forward-suppression factor sin 2 9 in (JHJ) 
is ignored. For high frequencies this has no effect. For low frequencies, the differences at 
small angles only plays a role for showers nearly parallel to the surface normal, while the 
effects of changing the normalization near the Cerenkov angle is important also for more 



horizontal showers. A measure of the effective angular spread Ac of the emission around 
the Cerenkov angle 9 c is given in terms of 



ln (2) 1 ™ 
For the example in Fig. [2] this becomes Ac — 44°. The Cerenkov radio waves are expected 



to be 100% linearly polarized [2JJ. Radio signals should therefore be observed with both the 
LCP and RCP modes. A simultaneous triggering of both the LCP and RCP modes (with 
50% of the total intensity in each mode) will be a good signature of the transient Cerenkov 
radio wave emission from UHE particles. In addition, the large number of GMRT dishes 
should allow pointing discrimination and spatial-temporal coincidence measurements to be 
used for signal identification and background rejection. 

III. APERTURE FOR COSMIC RAY AND NEUTRINO EVENTS AT GMRT 

For cosmic rays, which produce showers by hadronic processes, the shower energy E a is 
equal to the energy of the primary particle. When instead the showers are produced by 
deep inelastic neutrino scattering off target nucleons, the fraction of the total energy which 



produces the hadronic showers is, on average, E s = 0.25Ep 23]. 

We model the lunar regolith as a homogeneous medium with dielectric permittivity e = 3, 
hence refractive index n = 1.73. The Cerenkov angle is 9c = 54.7°. The radio waves are 
also absorbed in the regolith, which we take to have a frequency dependent attenuation 
length A r = 9/(z//GHz) m |22j. This means the attenuation length of radio waves is always 
much smaller than the mean free path of the neutrinos, conveniently parameterized as X u = 
670(1 EeV / E u ) ' 363 km using the cross sections of [23[]. As can be seen from this expression, 
the mean free path for an UHE neutrino is not sufficient to escape the Moon. Most of the 
passing neutrinos will therefore interact in the material. 

The threshold energy required for showers to be detectable is determined by the sensitivity 
of the radio telescope. This is commonly presented in terms of a noise intensity 

N VSt^A^ 

where T sys is the system temperature, ks is Boltzmann's constant, 5t the integration time, 
Au the bandwidth, and A e s the effective area of the telescope. In units of Jansky (1 Jy = 



10 26 W/m 2 /Hz) the noise intensity is given by 

F N = 2.76 x io 3 _J5^ . (8) 

v / 5tA^(A eff /m2) 

At GMRT, for frequencies 150, 235 and 327 MHz we use the effective area A e ^ = 30000 
m 2 while for the higher frequencies A c q = 18000 m 2 . The integration time 8t = 1/A U for 
the bandwidth limited case. In Table 1 we list the system temperatures at the different 
observation frequencies and the corresponding noise levels. Using equation we can solve 
for E s at the threshold required for measurement with the radio telescope (obtained for 
9 = 9c and F = -F/v). If we take a required signal-to-noise ratio er, the threshold shower 
energies E^ which can be measured at GMRT at the different observation frequencies are 
listed in Table 1. 

The radio waves produced below the lunar surface get attenuated while propagating to the 
surface and they get refracted at the surface before reaching the telescope. The transmission 
coefficients are different for the electric field polarization which is parallel and perpendicular 
to the lunar surface. It is therefore convenient to express the Cerenkov radiation intensity 
(J1J) in terms of the electric field magnitude £ of the radio waves. The conversion is 

F 10 7 / 8 \ 2 Av 



Jy 24vr V/uV/m/MHz J MHz ^ 
Using this conversion we also list the threshold electric fields measurable at GMRT in Table 
1. If attenuation is factored in separately, the quantity RE remains a constant between the 
point of emission and the point of refraction at the surface: 

-zp ( E s \ w ( ^ \ A*V 



Jg = - 2 ^ ^" UiW^ fc (10) 

At the surface, there is partial reflection which changes the intensity of the outgoing electric 
field. The electric field observed at the telescope is determined relates to the incident electric 
fields by the Fresnel relations 

R£j_ 2 cos r 



r, 



(R£±)mc n cos % + cos r 



21,= ,^' = 2C ° Sr (11) 
(i?£||)inc ncosr + cosi 

where the angles of incidence i and refraction r are related by Snell's law. The refraction 
angle is fixed by the polar angle 9 S through the geometrical condition 

Mr) = m • (12) 



Here R(9 S ) = [R 2 M + (Rm + h) 2 — 2(R M + h) R M cos 9 S ] 1 ^ 2 is the distance between the detector 
and the point of refraction on the lunar surface, and h = 384 000 km is the average Earth- 
Moon distance. Putting everything together, the parallel component of the electric field at 
the telescope on Earth is given by 

£|| = ^W)0R£) inc , (13) 

where (RS)i nc is given by ffTUj) as a function of 9 and E s . The dependence on the polar angle 
9 S of the radio flux on Earth is shown in Fig. [31 

Evaluating (_R£) inc at 9 = 9c and equating S\\ with the threshold field £ t h which can be 
measured with the telescope, we obtain the threshold energy E t h as a function of the polar 
angle 9 S . Fig. H] illustrates that for observing the radio waves from the limb of the Moon 
(9 S — > 7r/2 for h ^> Rm)-, one has to have a higher threshold energy of the shower owing 
to the fact that the transmission coefficients T± and Tji vanish when r — > n/2. When h is 
not large compared to Rm (which is the case for a telescope aboard a satellite orbiting the 
Moon) then only a small disc around the center of the Moon will be visible. 

By equating the signal £u = yfaE^ at the Cerenkov angle we determine the threshold 
Eth{9 s = 0) at different frequencies and bandwidths. The values so obtained are listed in 
TableE 



v (MHz) 


T 


(K) 


Ai/ (MHz) 


f n (Jy) 


£ N (^V/m/MHz) 


£th(0)/v^ (eV) 


150 


482 




40 


44 


.0029 


1.28 x 10 20 


235 


177 




40 


16 


.0017 


5.04 x 10 19 


325 


108 




40 


10 


.0013 


2.9 x 10 19 


610 


101 




60 


15 


.0014 


1.68 x 10 19 


1390 


72 




120 


11 


.0008 


5.6 x 10 18 



TABLE I: GMRT parameters, sensitivity and threshold energy. Fn is the noise intensity of GMRT 
and £_\f the corresponding electric field. The threshold energy in the last column is calculated with 
a = 25. 



The event rate that would be expected at the telescope can be related to an isotropic 
flux $ of UHE particles on the Moon through 



where i = {CR, u} denote the type of primary particle and Ai(E) is an aperture function 
corresponding to the effective detector area. We will differentiate between neutrinos, which 
have a large mean free path and therefore can penetrate the interior of the Moon before 
producing the Cerenkov radio waves, and the strongly interacting cosmic rays (referred to 
by the subscript CR) which can penetrate only small distances below the surface of the 
Moon. The aperture can be further decomposed into an angular aperture AQi(E, 9 S ) and a 
geometric area factor for the Moon 

Ai(E) = 2irR 2 M J Atti(E,0 s )dcos6 s (15) 

with Rm = 1760 km. To evaluate the aperture, we use the analytical methods described in 
For the case of cosmic rays, the angular aperture is given by 

AQ C r(E,9 s ) = J cos/3e[S(E,e s )-S th ] x 6(cos/5)dadcos/3, (16) 

where (3 and a are the polar and azimuthal coordinates of the ray normal to the Moon's 
surface in a system where the shower direction defines the z axis. The full geometry and the 
different angles are described in Fig. [U Using the Gaussian approximation for the radiation 
pattern, the angular aperture for cosmic rays is obtained by directly integrating over the 
shower coordinates (9, 6): 



/■sin(i) M m 

An CR = 2 dcos6» / d0[sin(0) sin(i) cos(0) - cos(0) cos(i)]. (17) 

Jl/n-D JO 



Here cos(0 m ) = cos W cos W anc i j n the lower integration limit D = w-ln-jP-. This takes 

sin(i) sm(6) ' y k i?th 

into account all directions with sufficiently strong radio emission. After performing this 
integration, the expression for the angular aperture reduces to 

_j / cos 9 V 



Att C R = cos — ; — — I - sin( 
sin 2 



sm(9 max ) cos(«)0 max + cos(^ max ) smM sm(0 max 



(18) 

where cos(0 max ) = 1/n — D and cos(0 max ) = sin(e max j • ^ ne aperture for cosmic rays 
is obtained by substituting ( fl8l) in ( fl5l) and integrating over the polar angle 9 S . 

When the UHE primary is instead a neutrino, it can produce showers deep below the 
surface of the Moon and there will be considerable attenuation of the radio waves which 
travel distances longer than A r below the surface. For the neutrino induced showers, the 



S 



aperture is denned in the same way as for the CR, but the angular aperture is now given 



m by 

AQ U (E,6 S ) = J y J <S>{£(E,9 S ,6) exp[-z/(2X r cosi)]-£ th }xe^p[-L(z,(3)/X v ]xdadcos(3, 

(19) 

where L(z, (3) is the distance the neutrino travels inside the material to reach the interaction 
point at a distance z below the surface. In performing this integration we allow z to go below 
the known depth of the regolith. The aperture can therefore pick up contributions from 
sufficiently strong signals coming from deep showers, especially for the lower frequencies. 
Numerically we find for the worst case (when v = 150 MHz), that imposing a sharp cutoff 
at a depth of 20 m would reduce the aperture by nearly an order of magnitude, similarly to 
what was discussed in [17 {. 

It has been observed [18( , that if the absorption length of radio waves A r is much smaller 
than the neutrino mean free path X u (which is indeed here the case) then Eq. (fl9l) reduces 
approximately to 

AQ U (E,9 S ) = j ln(S(E)/£ th )e[£(E,6 s ,9) - £ th ) x dadcos/3. (20) 



As for the cosmic rays, the total aperture is obtained by substituting (fl9j) into (1151) and 
integrating over the polar angle 9 S . 

To estimate the sensitivity of GMRT to cosmic ray and neutrino events we have evaluated 
the angular apertures by employing this technique and performing numerical integrations 
for the different parameters given in Table [B In the next section we will discuss these results 
further in the context of prospective flux limits. 

IV. LIMITS ON THE FLUX OF COSMIC RAYS AND NEUTRINOS 

Should no events be observed at GMRT during observation over a time T, an upper limit 
can be established on sufficiently smooth UHECR and neutrino fluxes at the Moon. The 
conventional model-independent limit |26J] is given by 

2 d$i Ei 
E *&E v - Su »A l {E s = y l E l )r (21) 

where still i = {z/, CR}, yen = 1 and y v = 0.25. The Poisson factor s up = 2.3 for a 
limit at 90% confidence level. The limits on the flux of UHECR that could be established 



9 



for 100 hours and 30 days of observation time at GMRT are shown in Fig. [5] and Fig. [6] 
respectively. We also show the results of carrying through our calculations for the LOFAR 
parameters given in [l?j]. When compared to the LOFAR Monte Carlo simulation results, 
the two methods of calculation agree within a factor two over the full energy range. This 
is an acceptable discrepancy in level with known uncertainties from e.g. the regolith depth. 
A residual difference of ~ 30% is also expected since we used the Gaussian approximation. 
The limit we obtain for the LOFAR parameters is less stringent than the one published, so 
in this respect our result constitutes a conservative estimate. ^From comparing the limits 
for different frequencies, it can be seen that low frequency observations give more stringent 
limits on the flux at the expense of a higher threshold. This is due to the well-known increase 

m 

in the aperture 11711 from radiation spreading at lower frequencies. The figures also show 
the Auger data [2J] on UHE cosmic rays. Although the extrapolation to higher energies is 
highly uncertain, the GMRT would most probably be sensitive to a post-GZK proton flux. 

Similarly for the UHE neutrinos, prospective limits on their flux for T = 100 hours 
and T = 30 days are shown in Figs. H and El Also here we show a calculation for the 
LOFAR parameters, again in quantitative agreement with previous results. Since many 
radio experiments exist for UHE neutrino detection, we have compiled a comparison in 
Fig. [9j This figure contains, in addition to the GMRT results for v = 150 MHz with two 
times, the already existing limits from RICE [25J, GLUE (lj, FORTE 



different observation 
3 and ANITA-lite [7 
calculated for ANITA 



. Also we have indicated t he p rospective future limits that has been 



7], LOFAR (l7| or LORD 



V. SUMMARY AND CONCLUSIONS 



We have calculated the potential for GMRT to detect UHE cosmic rays and neutrinos 
through their radio wave emission produced when showering in the lunar regolith. Our 
results indicate that GMRT could be competitive to future experiments in the E > 10 20 
eV range. If one assumes that the CR spectrum continues beyond the GZK limit with 
unchanged energy dependence, observation of these particles with the GMRT should indeed 
be possible. For UHE neutrinos there exists a theoretical upper bound on the flux from 
cosmogenic sources, as given by Waxman and Bahcall [27J: < 2 x 10 8 GeVcm 2 

g -i sr -i The GMRT could, using a mere 30 days of observation time, probe fluxes a factor 



10 



five smaller. As a benchmark scenario, we indicate in Fig. [9] the predicted flux from one TD 



model 28J where the mass scale resides around 10 22 eV. 

It is notable that the GMRT (for the low frequencies) has only a somewhat higher thresh- 
old, and slightly worse sensitivity for this type of experiment, than will the LOFAR tele- 
scope. This clearly points to the interesting potential of utilizing the GMRT for UHE 
particle searches in the near future. The results we have presented here obviously depend 
on to what extent the experimental realization of these measurements are possible at the 
GMRT facility. We therefore foresee a future analysis taking more thoroughly into account 
the requirements on the technical infrastructure and the experimental techniques for signal 
identification and background discrimination. 
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FIG. 1: Geometry of an UHE CR or neutrino event which generates Gerenkov radiation of radio 
waves in the lunar regolith. The Earth is located at a mean distance h from the lunar surface. 
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FIG. 2: Radio flux at Earth from the center of the Moon vs. emission angle at v = 150 MHz 
and bandwidth Av = 40 MHz for a shower initiated with E s = 10 20 eV. The dashed line contains 
the complete ^-dependence in Eq. ([5]), whereas the full line shows the result of using the Gaussian 
approximation . 
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FIG. 3: Radio flux at Earth vs. the lunar polar angle 6 S for 6 = 6 C at v = 150 MHz and bandwidth 



Au = 40 MHz for a shower initiated with E s = 10 20 eV. 
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FIG. 4: Shower threshold energy vs. polar angle 9 S on the Moon for a = 25. 
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FIG. 5: Model independent limits on UHECR flux at different 
vation time with GMRT. Auger data points reproduced from 
comparison. 
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FIG. 6: Limits on UHECR flux at different frequencies for 30 days of observation time with GMRT. 
Auger data as before. 
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FIG. 7: Model independent limits on UHE neutrino flux at different frequencies for 100 hours 
of observation time with GMRT. For comparison we show also a limit calculated for the same 
observation time with the LOFAR parameters of [it| . 
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FIG. 8: Limits on UHE neutrino flux at different frequencies for 30 days of observation time with 
GMRT. LOFAR calculation with the same parameters as before but also for T = 30 days. 
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FIG. 9: Prospective flux limits on UHE neutrinos from GMRT shown for effective exposure times 
of 100 hours and 30 days. The current best limits from radio experiments ANITA-lite [7j, GLUE 



12]. FORTE 3. aild RT( ' E 

experiments ANITA [7], LOFAR 



25]are shown. For comparison the expected limits from future 
17] and LORD []j] are also included. The WB line indicates the 



theoretical upper limit of Waxman-Bahcall 



271 ] on the cosmogenic neutrino flux. TD refers to the 



Topologi cal Defect model with Mx = 2 x 10 22 eV described in 



28(]. 
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